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I report on some recent developments in nonperturbative studies of finite temperature 
QCD with dynamical quarks on the lattice. I discuss new studies of improved lattice actions 
and their application to finite temperature QCD. I also summarize the status of lattice 
investigations about the order of the finite temperature QCD transition for the case of two 
flavors of degenerate light quarks, using both staggered and Wilson lattice fermions. 

o' 

£f~) • §1. Introduction 

(N 

The finite temperature transition in QCD is characterized by recovery of the 

Q^ ■ spontaneously broken chiral symmetry and also by breakdown of the confinement at 

temperatures higher than a critical value T c . Because both of these characteristic 

Q_i! properties of QCD are essentially non-perturbative in nature, the most systematic 

O ■ way of studying the transition is to investigate it on the lattice. 

A large number of numerical simulations have been performed to clarify the 
nature of the QCD transition from the first principles of QCD. In order to make a 
precise prediction for the real world from the results obtained on finite lattices, we 
perform an extrapolation to the limit of vanishing lattice spacing, the continuum 
limit. For a reliable extrapolation, we have to simulate on sufficiently fine lattices. 
On the other hand, in a study of finite temperature QCD, the spatial lattice size 
should be sufficiently large in order to approximately realize the thermodynamic 
limit. With the limitation of the currently available computer power, this means 
that we have to suppress the lattice size in the time direction iVj; currently Nt ~ 4-6 
is used for major simulations of QCD with dynamical quarks (full QCD). Because 
the temperature on the lattice is given by T = l/iVja with a the lattice spacing, the 
corresponding lattices become rather coarse; for T ~ T c ~ 100-200MeV, a ~ 0.2- 
0.4fm. Subsequent lattice artifacts sometimes make the analysis and interpretation 
of lattice results not straightforward. 

Here, recent developments in improved lattice actions opened us a possibility 
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2 K. Kanaya 

to reduce the lattice artifacts already on coarse lattices and, therefore, to extract 
accurate results for finite temperature QCD transition with the present computer 
power. In section g, I report these developments in improved lattice actions and 
their application to finite temperature QCD. 

In a phenomenological study of quark-gluon plasma in heavy ion collisions and 
in the early Universe, the most fundamental information required is the order of the 
QCD phase transition. Recently, there have been several developments in lattice 
QCD on this issue, especially for the case of two flavors of degenerate light quarks. 
In sect. [|, I report about the status of these recent studies using both staggered and 
Wilson lattice fermions. This also includes one of the first fruitful applications of 
improved actions to the study of QCD at finite temperatures. 

For other developments in lattice QCD, I refer reviews in recent Lattice confer- 



ences HI 

§2. Improvement of the lattice action 

An action on the lattice is chosen such that the continuum action is recovered 
when we let the lattice spacing a — ► in the action. For the gluon part of the 
action, the minimal choice is the standard one-plaquette action by WilsoncP, which 
is the most local and geometrically the simplest lattice gauge action. On the lattice, 
however, we have freedom to introduce less local terms to the action without affecting 
the continuum limit. Here, the speed of approach to the continuum limit does depend 
on the choice of the action. Therefore, a judicious choice of the action can suppress 
lattice artifacts in physical observables even at moderate values of the lattice spacing. 
Such actions are called "improved actions". Importance of improvement is much 
stressed recently (for reviews see Refs.EP'l§). 

The basic idea behind improvement may be obtained by considering the lattice 
derivative; the naive derivative A x f(x) = £:[f(x + a) — f(x — a)] converges to 
the continuum derivative with an error of 0(a 2 ). We can reduce this error down 
to 0(a 4 ) by replacing A x — ► A x — \A X which operates on fields at up to next 
nearest neighbor points. Similar substitution is effective also to obtain a lattice 
action that approaches to the continuum action much smoothly. The life is, however, 
not so easy, because what we want to obtain is not a smoother lattice action, but 
an action which leads to smaller discretization errors in the physical observables 
(P) = Z^ 1 J T><pP[<p]e~ Wi which can be obtained only after a non-perturbative 
calculation, in general. Many different strategies have been proposed to obtain such 
S. Two major approaches are the renormalization group methods proposed first 
by Wilson, and the perturbative improvements by Symanzik. In any case, improved 
lattice actions contain additional terms (interactions that usually have a wider spatial 
size) compared with the standard action. Because such additional terms makes the 
simulation quickly difficult and time-consuming, the efficiency of improvement should 
be tested for each of the new terms introduced. 

2.1. Comparison of improved actions 
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Fig. 1. Static quark potential obtained by a full QCD simulation performed on 12 x 32 lattices 
using various lattice actions. EP The coupling parameters /3 — 6/g 2 and K are adjusted such that 
we obtain the lattice spacing a ~ 0.2fm and the ratio of the pseudo-scalar and vector meson 
masses mps/my ~ 0.8. 



Improved lattice actions have been studied mainly in quenched QCD (QCD 
without dynamical quark loops). In finite temperature QCD, however, dynamical 
quarks are essential. Quite recently, the CP-PACS Collaboration have performed a 
systematic investigation of how various terms added to the gauge and quark actions, 
taken separately, affect light hadron observables in full QCDEP. A new dedicated 
parallel computer CP-PACS cP is used for the simulation. The effect of improvement 
is studied comparing, for the gauge part of the action, the standard one plaquette 
action (P) and an improved action (R) that contain 1x2 rectangular loops with 
the coefficients obtained by a renormalization group studylZP, and for the quark 
action, the standard Wilson quark action (W) and an improved Wilson quark action 
containing the so called clover termtj (the clover action: C). 

Simulations are carried out in two-flavor QCD. Because the extent of improve- 
ment to be clearer at a coarser lattice spacing, we adjust the gauge coupling param- 
eter 6/g 2 such that the lattice spacing becomes a ~ 0.2fm. 

Figure 1 summarizes the results for the static quark potential using the four 
different combinations for the action (PW, PC, RW, and RC). Different symbols 
correspond to potential data measured in different spatial directions on the hypercu- 
bic lattice, along the vector indicated in the figure. The jagged result for the standard 
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Fig. 2. Octet baryon mass as a function of PS meson mass squared, normalized by vector meson 
mass at a ~ 0.2fm using various lattice actions. W The diamond near the lower left corner is the 
experimental point computed with N, n and p. 



action (PW) means that the rotational symmetry is sizably violated on this coarse 
lattice. A drastic improvement of the rotational symmetry is seen when we replace 
the gauge action with the R action, as observed previously also in the quenched 
casecP. In this regard, the quark action has much less effect. Fig. 2 shows the results 
for light hadron masses. The solid curve represents the result of a phenomenological 
mass formula by OnolI2P , which is approximately reproduced also on the lattice when 
a is sufficiently small (a < O.lfm for the PW action). We see that improvement of the 
quark action makes the spectrum quite close to the phenomenological curve already 
at a « 0.2fm. 

These results are very encouraging: They show that, even with the presently 
available computing power, it may be possible to obtain accurate results near the 
continuum limit also in full QCD, when we use an improved action. 

2.2. Improved actions for full QCD at finite temperatures 

One of the first systematic applications of improved action to the study of finite 
temperature QCD has been done by the Tsukuba group Ef' HIP. Their motivation to 
adopt an improved action is to remove severe lattice artifacts encountered in a finite 
temperature simulation of QCD with the standard Wilson quark action EJ)>E-§. Using 
the RW action discussed in the previous subsection, these lattice artifacts are shown 
to be well removed. Further physical outputs from this study will be discussed in 
Sec. 3.2. 

Through this study, it became clear that improvement of the lattice action is 
essential in order to study several important topics in finite temperature QCDH3. 
Other major research groups have also began to study the finite temperature QCD 
with improved actions li^P . 

The MILC Collaboration combined a Symanzik type improved gauge action and 
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the clover quark action C, and found good improvement liZP similar to that observed 
in the case^of the RW action. Comparing these results together with an unpublished 
PC dataEJ), where the lattice artifacts are found to be remaining, we conclude that 
the big reduction of the lattice artifacts noted in the simulation of finite temperature 
QCD is mainly due to improvement of the gauge action. The MILC Collaboration 
studied also an improved staggered quark action, aiming at a better flavor symmetry 
on a coarse lattice. 

The Bielefeld group also extended their study of improved pure gauge theories I12P 
to include improved staggered and Wilson quark actionsS . One of their motivations 
is to achieve a faster convergence of thermodynamic quantities to the large Nt limit: 
Thermodynamic quantities have errors of 0{1/Nf) on the lattice because the relevant 
fluctuations include those with spatial size ~ 1/T = N-ta, that feel discretization 
errors when Nt is not large enough. These errors are automatically removed when 
we take the continuum limit a — > keeping the temperature T = 1/N t a fixed, i.e., 
Nt — > oo simultaneously. [Practically, N t >8 is required to obtain thermodynamic 
quantities at T ~ 0(T C ) to an accuracy of several percent in quenched QCD using 
the standard one-plaquette action.] Because an improved action has improved short- 
distance properties, we expect these errors to be also smaller with improved actionsJiil 

In spin models and also in quenched QCD, good improvements are reported using 
a class of improved actions called fixed point actions ("perfect actions"), which are 
designed to be at the UV fixed point of a renormalization group transformation and 
therefore to remove all lattice artifacts when infinite number of coupling parameters 
3dfP. Trials to construct a fixed point action for full QCD are also 
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§3. Order of the chiral transition in two-flavor QCD 

Let us now turn our attention to the topics of the order of of the QCD transition 
from the high temperature quark-gluon-plasma phase to the low temperature hadron 
phase, which is supposed to occur at the early stage of the Universe and possibly at 
heavy ion collisions. It is, in particular, crucial to know whether the transition is a 
first order phase transition or a smooth transition (second order phase transition or 
crossover). 

The nature of the QCD transition is considered to depend on the number of 
quark flavors Np sensitively. The physically interesting cases are the case of two and 
three flavors of degenerate light quarks (Np = 2 and 3), and also a more realistic 



*' The Bielefeld group extended the criterion to include improvement also in the "high tem- 
perature limit" on lattices with a fixed small Nt ~ 4: Since a is a decreasing function of 6/g 2 , a 
larger 6/g 2 corresponds to a larger T = 1/Nta. For fixed Nt, however, the "high temperature limit" 
6/g 2 —* oo corresponds to a temperature higher than T = oo in the usual procedure in which the 
limit T — » co is taken only after the limit a — » with fixed T. Because this includes an exchange 
of the limiting procedures, 1 think this extended criterion to be slightly different in nature from the 
usual criterions for improvement concerning the convergence to the continuum limit. The difference 
becomes important when the number of adjustable coupling parameters in the improved action is 
restricted. 
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case of two light and one heavy quarks (Np = 2 + 1). The case Np = 2 corresponds 
to the case where the third quark s is much heavier than the relevant energy scale 
for thermal processes near the critical temperature, m s 3> T c , while the case Np > 3 
corresponds to the case in s <C T c . 

Using universality hypothesis, the nature of the finite temperature QCD tran- 
sition near the chiral limit can be studied by a Ginzburg-Landau effective theory 
respecting the chiral symmetry of QCD, the effective a modeled. For Np > 3, a 
first order transition is predicted. For Np = 2, on the other hand, the order of 
the transition is not quite definite in the effective a model; a first order transition 
is predicted when the anomalous axial Ua(1) symmetry is effectively restored at 
the transition temperature, while a second order transition is expected otherwise. 
Because the Ua(1) breaking operator is a relevant operator whose coefficient grows 
towards the IR limit under a renormalization group transformation, the transition is 
more likely to be second order EJ. A non-perturbative study is required to determine 
the order of the transition conclusively. 

Therefore, understanding the nature of the QCD transition for Np = 2 is an 
important step toward the clarification of the transition in the real world. 

When the chiral transition is second order, we expect that the transition turns 
into an analytic crossover at non-zero m q , while when the chiral transition is first 
order, it will remain to be first order for small m q . However, it is difficult to nu- 
merically distinguish between a very weak first order transition and a crossover, 
especially at small m q : In order to confirm the expected crossover numerically, we 
have to study the lattice size dependence to see if the formation of singularity (e.g. 
the increase of the peak height of a susceptibility with increasing the lattice volume) 
stops on sufficiently large lattices. This becomes more and more difficult at small 
m q . 

Here, the universality provides us with useful scaling relations that can be con- 
fronted with numerical results of QCD, in order to test the nature of the transition: 
It is plausible from an effective a model that, when the chiral transition is of second 
order, QCD with two flavors belongs to the same universality class as the three di- 
mensional 0(4) Heisenberg modelEP. The 0(4) model is much simpler than the a 
model, and its scaling properties are well studied. For example, at small external field 
h near the critical temperature T c for h = 0, the pseudo-critical temperature T pc (h) 
and the peak height of the magnetic and thermal susceptibilies follow T pc — T c ~ h Zg , 
X™ ax ~ h~ Zm , and X r ax ~ h ' Zt , wh ere z g = 1/(35, z m = l- 1/5, and z t = (1 -/3)//36 
in terms of the 0(4) critical exponents /3 and 5. Here the values of j3 and 5 for the 



0(4) model are well established S*. In QCD, we identify T ~ 6/g 2 , h ~ m q , and 
M ~ (&&). 

In lattice QCD, an additional complication should be noted because no known 
lattice fermions have the full chiral symmetry on finite latticescB'. Two conven- 
tional lattice fermions are the staggered and Wilson fermions. In the formulation 
of staggered fermionsEIP, the flavor-chiral symmetry S\J(Np)i x S\J(Np)r x U(l) 
of massless iV^-flavor QCD is explicitly broken down to \J(Np/4)i x \J(Np/4:)ji at 
a > 0. Moreover, this action is local only when Np is a multiple of 4. For the 
physically interesting cases Np = 2 and 3, a usual trick is to modify by hand the 
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Fig. 3. (a) Peak height of the magnetic susceptibility for Nf — 2 QCD with staggered quarks as a 
function of the spatial lattice volume iVf . (b) The same data as a function of the lattice volume 
rescaled by zero-temperature pion correlation length. 



power of the fermionic determinant in the numerical path-integration. This neces- 
sarily makes the action non-local, that poses conceptually and technically difficult 
problems. With the Wilson fermion actionS, on the other hand, the flavor symme- 
try is manifest also on the lattice, and the action is local for any Np. This feature is 
quite attractive in a study of the Np dependence of the QCD transition. However, 
the chiral symmetry is explicitly broken at a > 0, that requires additional tedious 
analyses of chiral properties. 

These broken symmetries of lattice fermions are expected to restore in the con- 
tinuum limit. On a coarse lattice used in a finite temperature simulation, however, 
we may encounter sizable deviations from the scaling behavior expected in the contin- 
uum limit. Therefore, the appearance of the 0(4) scaling is also a useful touchstone 
to test the recovery of the chiral symmetry on the lattice when the chiral transition 
is of second order. 

3.1. Results with staggered quarks 

The 0(4) scaling was first tested on the lattice for staggered quarks by the 
Bielefeld groupEH*. Based on simulations on an 8 3 x 4 lattice at m q a = 0.02, 0.0375, 
and 0.075 using the standard action, they obtained z* = 0.77(14), z m = 0.79(4), 
and z t = 0.65(7), where the corresponding 0(4) values") are 0.537(7), 0.794(1), and 
0.331(7). The result for z m is consistent with the 0(4) value while other exponents 
are in disagreement with the 0(4) values. 

Possible causes of the discrepancy are (i) m q is not small enough to see the criti- 
cal behavior in the chiral limit, and (ii) the spatial lattice volume is not large enough 
to obtain the observables in the thermodynamic limit. Two additional caveats are in 
order for Np = 2 staggered quarks: (iii) The symmetry in the chiral limit at a > 
is 0(2) instead of 0(4). Practically, however, the values of the 0(2) exponents are 
almost indistinguishable from the 0(4) values with the present numerical accuracy, 
(iv) The action is not local. Therefore, an assumption behind the universality ar- 
gument can be violated so that some non-universal behavior may appear 113. The 
correct continuum chiral limit with the 0(4) symmetry will be obtained only when 



8 K. Kanaya 

we first take the continuum limit a — > and then take the chiral limit. In addition 
to these points, we also have to check technical details in the numerical simulation; 
the accuracy of the methods to simulate the system, such as the the finite step-size 
error and the dependence on the convergence criterion for fermion matrix inversion. 

A systematic study of the quark mass dependence as well as the lattice volume 
dependence is in order. We performed a series of simulations on 8 3 x 4, 12 3 x 4, and 
16 3 x 4 lattices at m q a = 0.01, 0.02, 0.0375, and 0.0750. The Bielefeld group also 
extended their study to larger spatial lattices EiP . The results obtained are consistent 
with each other. It turned out that determination of critical exponents on 8 3 x 4 
lattices suffers from a sizable finite lattice-size effect for m q a < 0.0375. 

From the lattice-size dependence of the magnetic susceptibility Xm for fixed value 
of quark mass, shown in Fig. 3(a), we conclude that the transition is a crossover for 
m q a > 0.02; the peak height Xm ax f° r rn <? a = 0-02 stabilizes on spatial lattices larger 
than 12 3 . For m q a = 0.01, on the other hand, x™ ax 1S increasing up to our largest 
spatial lattice of 16 3 . If this increase is maintained up to infinite volume, then we 
have a first order transition at this quark mass. However, examining the lattice 
volume dependence of Monte Carlo time histories and histograms, we could not find 
any clear indication of a first-order transition at m q a = 0.01. Furthermore, when 
we rescale the lattice volume by zero-temperature pion correlation length, we find 
that the lattice volume 16 3 for m q a = 0.01 approximately corresponds to the volume 
12 3 for m q a = 0.02, where the increase of x™ ax terminates [Fig. 3(b)]. Therefore, it 
is possible that the increase of Xm ax fo r Trina = 0.01 seen in our data is a transient 
effect. 

Assuming that the finite size effect is sufficiently small on the 16 3 x 4 lattice, 
we fit the data at the four values of m q a. We find z g = 0.64(5), z m = 1.03(9), 
and Zt = 0.82(12). (Removing the data for m q a = 0.01 gives slightly smaller but 
consistent values with larger errors.) The results for z g and z m sizably deviate 
from the 0(2) or 0(4) values. We also studied other exponents related to energy 
expectation value, which also largely deviate from 0(2) and 0(4) values. On the 
other hand, the identity z g + z m — z t = 1 expected for a second-order fixed point with 
two relevant operators is approximately satisfied. Thus, the exponents are consistent 
with a second-order transition at m q = 0. Our study of the scaling function also 
shows a similar trend: With measured values of exponents, data for the susceptibility 
Xm exhibits a reasonable scaling as a function of j3 and m q a, although results are 
much worse when we adopt the 0(4) exponents. 

In summary, we find the determination of the nature of the two-flavor chiral 
transition with staggered quarks using the standard action to involve subtle prob- 
lems. While our data so far do not contradict a second-order transition at m q = 0, 
the exponents take quite unexpected values, at least in the range m q a > 0.01. Ev- 
idently further work, possibly on larger spatial sizes and smaller quark masses, is 
needed to clarify this important problem. 

3.2. Results with Wilson quarks 

We now study the issue using Wilson quarks. It turned out that Wilson quarks 
in the standard action lead to several unexpected phenomena on lattices with iVj = 4 
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Fig. 4. The Polyakcm. loop and the pion screening mass obtained with Wilson quarks with a RG 
improved actiontj'. In these figures, the conventional notation j3 = 6/g 2 is used, not to be 
confused with a critical exponent. Larger /3 corresponds to a higher temperature. The horizontal 
axis 1/K — 1/Kc is proportional to m q a. 



and 6: On these lattices, the transition becomes once very sharp when we increase 
m q from the chiral limitEj^O, in contrary to the expectation in the continuum limit 
that the chiral transition becomes weaker when we increase m q . Together with other 
strange behaviors of physical quantities nearthe transition point, this phenomenon 
is identified as an effect of lattice artifactsEiP . Therefore, we apply an improved 
action. 

With the RW action discussed in Sec. 2.1, we first-find that the lattice artifacts 
observed with the standard action are well removed 0)'EJ). We also find that the 
physical quantities are quite smooth around the transition point at m q > 0, as 
shown in Fig. 4. The straight line envelop of m 2 at finite temperature (Nt = 4) 
shown in Fig. 4(b) agrees with m 2 obtained at low temperature (Nt = 8), and 
corresponds to the PCAC relation m 2 oc m q expected in the low-temperature phase. 
The smoothness of the physical observables strongly suggests that the transition is 
a crossover at m q > 0. 

Concerning the nature of the transition in the chiral limit, we find that the tran- 
sition becomes monotonically weaker when we increase 6/g 2 (see Fig. 4). Because 
the transition point shifts to larger 6/g 2 at larger m q , increasing 6/g 2 corresponds 
to increasing m q for the transition. We also note that m 2 in the chiral limit mono- 
tonically decreases to zero as we decrease temperature from above towards the chiral 
transition point, m 2 at the transition temperature for finite m q also shows a simi- 
lar monotonic decrease when we decrease m q . These results suggest that the chiral 
transition is continuous. 

For a more decisive test about the nature of the transition, a scaling study is 
required. From the universality argument we expect that magnetization M near the 
second order transition point can be described by a single scaling function: 



M/h l/d = f(t/h L ^ d ), 



\/$8\ 



(3-1) 



where h is the external magnetic field and t = [T — T c ]/T c the reduced temperature. 
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When the QCD transition is of second order in the chiral limit, the chiral condensate 
should satisfy this scaling relation with 0(4) exponents 1/(35 = 0.537(7) and 1/5 = 
0.2061(9)c3) and the 0(4) scaling function f(x). Our results for M are shown_in 
Fig. 5(a). We make a fit of M to the scaling function obtained for an 0(4) modelEHP, 
by adjusting [3 c t and the scales for t and h, with the exponents fixed to the 0(4) 
values. Figure 5(b) shows our result with x 2 /df = 0.61. The scaling ansatz works 
remarkably well with the 0(4) exponents. Our recent study shows that the situation 
holds also when we include data in the t < region 113 . On the other hand, a change 
of the exponents quickly makes the fit worse: For example, when we use the MF 
exponents suggested by Kocic and Kogut as a possibility for two-flavor QCDE3\ the 
data no more falls on the MF scaling function. 

The success of this scaling test with the 0(4) exponents strongly suggests that 
the chiral transition is of second order in the continuum limit. It also indicates that 
the chiral violation due to the Wilson fermion action is sufficiently small with our 
improved action, for the values of m q and 6/g 2 studied here. 

§4. Summary 

A simulation of finite temperature QCD sometimes suffers from sizable lattice 
artifacts caused by the coarseness of the lattice. Due to a requirement of large spatial 
lattice size for the thermodynamic limit, and also due to physical requirements to 
study a wide range of the parameter space, it is hard to remove these lattice artifacts 
when we use the standard lattice actions. Recent developments in improved lattice 
actions opened us a possibility to perform accurate simulations of finite temperature 
QCD with the present power of computers. A comparative study of improved full 
QCD actions shows that an efficient reduction of lattice artifacts can be achieved 
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already on lattices with a ~ 0.2fm, as used in major finite temperature simulations. 

Among the topics of finite temperature QCD, intensive lattice studies have been 
made recently on the nature of the chiral transition in QCD with two degenerate 
light quarks. Understanding it is an important step toward the clarification of the 
transition in the real world including the non-degenerate s quark. Numerical simula- 
tions show that the transition for Np = 2 is analytic crossover down to a very small 
value of the quark mass. In order to further determine the nature of the transition 
in the chiral limit, we have to study the scaling property and compare it with the 
0(4) scaling expected theoretically for the case of a second order transition. 

Using the standard action for the staggered quark, we find the determination of 
the nature of the two-flavor chiral transition to involve subtle problems. While our 
data do not contradict a second-order transition at m q = 0, the critical exponents 
take quite unexpected values. It is not yet clear whether this strange behavior is 
caused by large sub-leading terms in the scaling relations, or by the non- universality 
caused by the non-locality of the two-flavor staggered quark action. If the latter is 
the case, application of an improved action may solve the problem. 

A study of the issue using the Wilson quark has been difficult due to severe 
lattice artifacts encountered when we use the standard action. We find that these 
lattice artifacts are removed when we use an improved gauge action combined with 
the standard Wilson quark action. Furthermore, we find that the chiral condensate 
remarkably follows the scaling relation with 0(4) exponents and 0(4) scaling func- 
tion. This strongly suggests that the chiral transition in two-flavor QCD is of second 
order. 

This second order transition for the case m s = oo should turn into a first order 
transition when we decrease m s , as observed for the degenerate Np = 3 case. Because 
m s ~ 150 - 200 MeV is just of the same order of magnitude as the expected values of 
T c , we have to fine-tune the value of m s in order to study the nature of the transition 
in the real world. Unfortunately, the lattice simulations performed so far for the 
Np = 3 and Np = 2 + 1 casesEJ'tU' do not have that accuracy yet because, using 
standard lattice actions, we have sizable systematic errors due to lattice artifacts. 
At the moment, both possibilities of a first order transition and an analytic crossover 
are remaining for the physical point. 

Encouraged with the success of the scaling study for two-flavor QCD using an 
improved action.we began test simulations for Np = 3 and 2 + 1 applying the same 
improved actionEr. Many lattice groups began to simulate finite temperature QCD 
using improved actions. It may, however, take several more years to re-accumulate 
the basic data in a wide range of the many-dimensional full QCD parameter space, 
that is required to clarify the nature of the QCD transition in the real world. 

This work is in part supported by the Grants-in-Aid of Ministry of Education 
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